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Abstract 

Consider random matrices A, of dimension m x (m + n), drawn from an ensemble 
with probability density f(trAA'), with f(x) a given appropriate function. Break A = 
(B,X) into an rn x to block B and the complementary m x n block X, and define 
the random matrix Z = B~ 1 X. We calculate the probability density function P(Z) 
of the random matrix Z and find that it is a universal function, independent of f{x). 
The universal probability distribution P{Z) is a spherically symmetric matrix-variate t- 
distribution. Universality of P(Z) is, essentially, a consequence of rotational invariance 
of the probability ensembles we study. As an application, we study the distribution of 
solutions of systems of linear equations with random coefficients, and extend a classic 
result due to Girko. 
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1 Introduction 



In this note we will address the issue of universality of the probability density function 
(p.d.f.) of the product B~ 1 X of real and complex random matrices. 
In order to motivate our discussion, before delving into random matrix theory, let us 
discuss a simpler problem. Thus, consider the random variables x and y drawn from the 
normal distribution 

G(x,y) = ^e-^ . (1.1) 

Define the random variable z = ~. Obviously, its p.d.f. is independent of the width a of 
G(x,y), and it is a straightforward exercise to show that 

i.e., the standard Cauchy distribution. 

A slightly more interesting generalization of (jl.lj) is to consider the family of joint 
probability density (j.p.d.) functions of the form 

G(x,y)=f(x 2 + y 2 ), (1.3) 

where f(u) is a given appropriate p.d.f., subjected to the normalization condition 

oo 

f(u)du = - . (1.4) 

7T 


A straightforward calculation of the p.d.f. of z = | leads again to (|1.2j) . Thus, the 
random variable z = ^ is distributed according to Q1.2j) , independently of the function 
f(u). In other words, (|1.2|) is a universal probability density function. 1 P(z) is universal, 
essentially, due to rotational invariance of (|1.3|h More generally, P{z) must be 
independent, of course, of any common scale of the distribution functions of x and y. 
We will now show that an analog of this universal behavior exists in random matrix 
theory. Our interest in this problem stems from the recent application of random matrix 
theory made in to calculate the complexity of an analog computation process [2], 
which solves linear programming problems. 



2 The universal probability distribution of the product 
B X X of real random matrices 

Consider a real m x (m + n) random matrix A with entries 

Aia (i = 1, . . . m; a = 1, . . .m + n). We take the j.p.d. for the m(m + n) entries of A as 

G(A) = f(trAA T ) = f , (2-1) 



x We can generalize 11.31 somewhat further, by considering circularly asymmetric distributions G(x,y) — 
f{ax 2 + by 2 ) (with a,b > of course, and the r.h.s. of 11.411 changed to ^/aE/ir), rendering 11.21 a Cauchy 
distribution of width yjb/a, independently of the function f(u). 
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with f(u) a given appropriate p.d.f.. From 2 



J G(A) dA = 1 (2.2) 
we see that f(u) is subjected to the normalization condition 

oo 

/m(m + n) A 
u— l f(u)du = , (2.3) 

where 

is the surface area of the unit sphere embedded in d dimensions. This implies, in 
particular, that f(u) must decay faster than u ~ m ( m + n )/ 2 as u _> oo ; an d also, that if f(u) 
blows up as u — > 0+, its singularity must be weaker than u - m ( m + n )/ 2 . i n other words, 
/(u) must be subjected to the asymptotic behavior 

u m( m +n)/2 /(n) ^ q (2 _ 5) 

both as u — ► and it — > oo. 

We now choose m columns out of the m + n columns of A, and pack them into an m x m 
matrix B (with entries By). Similarly, we pack the remaining n columns of ^4 into an 
m x n matrix X (with entries Xj p ). This defines a partition 

A^(B,X) (2.6) 

of the columns of A. 

The index conventions throughout this paper are such that indices 

i, j, . . . range over 1, 2, . . . , m , 

p,q,... range over 1,2, ...,n, (2-7) 
and a ranges over 1, 2, . . . , m + n . 

In this notation we have tr AA T = Bfj + J2i, P X i P = ^BB T + trXX T , and thus (|2~Tj) 
reads 

G(B, X) = f(tiBB T + tiXX T ) . (2.8) 

We now define the random matrix Z = B~ 1 X. Our goal is to calculate the j.p.d. P(Z) for 
the mn entries of Z. P(Z) is clearly independent of the particular partitioning l|2.6|) of A, 
since G(B, X) is manifestly independent of that partitioning. The main result in this 
section is stated as follows: 



2 We use the ordinary Cartesian measure dA = d m ^ m+n) A = Y\ ia dA ia . Similarly, dB — d m B and 
dX = d mn X for the matrices B and X in and lITTHt . 



2 



Theorem 2.1 The j.p.d. for the mn entries of the real random matrix Z = B 1 X is 
independent of the function f(u) and is given by the universal function 

P{Z) = — ° T „ + „ , (2.9) 

[det(l + ZZ T )\ — 

where C is a normalization constant. 

Remark 2.1 The probability density function (|2.9f) is a special (spherically symmetric) 
case of the so-called 3 matrix variate i-distributions [HI E] : The m x n random matrix Z is 
said to have a matrix variate i-distribution with parameters M, E, £1 and q (a fact we 
denote by Z (g, M, E, Q)) if its p.d.f. is given by 

n m r / i i m\ n — i (m+n+q— 1) 

D(detE) _ 2(detn)-T [det (l m + E^^Z - M)^- X (Z - M) T jJ 2 , (2.10) 

where M, E and f2 are fixed real matrices of dimensions m x n, m x m and n x n, 
respectively. S and f2 are positive definite, and q > 0. The normalization coefficient is 

^=4. n ^ l 77^ . (2.H) 

It arises in the theory of matrix variate distributions as the p.d.f. of a random matrix 
which is the product of the inverse sqaure root of a certain Wishart-distributed matrix 
and a matrix taken from a normal distribution, and by shifting this product by M, as 
described in Our universal distribution H2.9j) corresponds to setting 

M = 0, E = l m , f2 = 1„ and q = 1 in (j2~TTl and ijOTTl . 

Remark 2.2 It would be interesting to distort the parent j.p.d. (|2,1|) into a non-isotropic 
distribution and see if the generic matrix variate t-distribution (|2.1U|) arises as the 
corresponding universal probability distribution function in this case. 

To prove Theorem (|2.1[) . we need 

Lemma 2.1 Given a function f(u), subjected to (|2.M|) . the integral 

I = JdBf{trBB T ) \detB\ n (2.12) 
converges, and is independent of the particular function f(u). 

Remark 2.3 A qualitative and simple argument, showing the convergence of (|2.12j) . is 
that the measure dfi(B) = dB \detB\ n scales as dfi(tB) = t m{m+n ^ d^{B), and thus has the 
same scaling property as dA in (|2.2|) . indicating that the integral ()2.12|) converges, in view 
of (|2.5|) . To see that I is independent of f{u) one has to work harder. 

3 Our notations in Remark 12,11 are slightly different from the notations used in In particular, we 
interchanged their E and f2, and also denoted their (T — M) T by Z — M here. Finally, we applied the 
identity det(l + AB) = det(l + BA) to arrive after all these interchanges from their equation (4.2.1) to 
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Proof. We would like first to integrate over the rotational degrees of freedom in dB. Any 
real m x m matrix B may be decomposed as [SJ |S] 



b = o x no 2 (2.13) 



rn ) ■ 



where 0\ 2 £ 0(m), the group ofmxm orthogonal matrices, and fi = Diag(wi, . . . ,u. 
where u>i, . . . ,uj m are the singular values of B. Under this decomposition we may write 
the measure dB as [31 |S] 

dB = dii(Oi)dfi(0 2 ) J] k 2 - oJ 2 j\d m uj , (2.14) 

i<j 

where dfj,(0\ 2) are Haar measures over the 0{m) group manifold. The measure dB is 
manifestly invariant under actions of the orthogonal group 0(m) 

dB = d(BO) = d{0'B) , 0,0'eO(m), (2.15) 

as should have been expected to begin with. 

Remark 2.4 Note that the decomposition (j2.13j) is not unique, since 0\D and T>0 2 , 
with T> being any of the 2 m diagonal matrices Diag (±1, • • • , ±1), is an equally good pair 
of orthogonal matrices to be used in (|2.13|) . Thus, as 0\ and O2 sweep independently over 
the group 0(m), the measure (|2.14[) over counts B matrices. This problem can be easily 
rectified by appropriately normalizing the volume V m = J d^{0\)d^{02)- One can show 4 
that the correct normalization of the volume is 

m(m+l) 

V m = ^ \ / ■ \ • (2.16) 

2 m nr=ir(i + |)r(|) 

Let us now turn to ()2.12|) . The integrals over the orthogonal group in (|2.12j) clearly factor 
out, and we obtain 

i=v m [ nh 2 -^i n<* n5>*)- (2 - l7) 

Joo i=l j<k \i=l ) \i=l ) 

Finally, we change the integration variables in (|2.17|) to the polar coordinates associated 
with the uJi. The angular part of that integral is fixed only by dimensionality and by the 
factor IIka; \ u j ~~ ^kl (11^1 ^i) n , and is thus independent of the function f{u). 
To prove that / < 00 we need only consider integration over the radius r 2 = J2iLi^f, 
since integration over the angles obviously produces a finite result. Using (J2.3|) . we find 
that the radial integral in question is 

00 00 

2 1 -1 o /* m(m+n) -. 



2 J *S'm(m+n) 



4 One simple way to establish 12.16L is to calculate J dB exp —^trB T B = (2n) 2 
V m j d m ui Yii<j ~~ w f I ex P — I X/i w i ■ r ^ ne ^ as ^ integral is a known Selberg type integral 0. 
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independently of f(u). I 

We are ready now to prove Theorem (|2.1j) : 
Proof. By definition 5 , 

P(Z) = J dBdX f{tiBB T + trXX T )5(Z - B^X) 

= J dBdX f(tvBB T + tiXX T )\detB\ n 5(X - BZ). (2.18) 

Integration over X gives: 

P(Z) = J dB f(tvBB T + trBZZ T B T )\detB\ n 

= JdBf[tvB(l + ZZ T )B T ] |detJ3| n (2.19) 

The m x m symmetric matrix 11 + ZZ T can be diagonalized as l + ZZ T = OAO T , where 
O is an orthogonal matrix, and A = Diag(Ai, . . . , A m ) is the corresponding diagonal form. 
Obviously, all Aj > 1, since ZZ T is positive definite. Substituting this diagonal form into 
(|2TT9~|t we obtain 

P(Z)= J dB f{tvBOAO T B T )\detB\ n . 

From the invariance of the determinant |det-B0| = |det-B| and of the volume element 
d(BO) = dB under orthogonal transformations we have: 

P(Z)=JdBf(trBAB T )\detB\ n . (2.20) 
Let us now rescale B as B = B\fX. Thus, 



det B = VdetA det B , and dB = (det A) f dB . (2.21) 
Finally, substituting (|2~2"T|) in (f2~2H)) we obtain 



dB , ~ ~ T\ ( IdetBl \" 
P(Z) = - —^ f(.tvBB T ) t^-^t 

(det A) 2 \ (det A) 2 / 

C C 



(2.22) 



(det A)( m +™)/ 2 [det(l + ZZ T )]( m + n )/ 2 ' 
where C is the normalization constant 

C = J dBf(trBB T )\detB\ n (2.23) 

rendering 

J P(Z)dZ = 1 . (2.24) 

C is nothing but the integral (|2.12jl . Thus, according to Lemma (|2.1[) . C < oo and is also 
independent of the function f(u). I 



5 Our notation is such that 5(X) = f]™ 1 ]J n =1 S(X lp ) = J]p=i (5(m) (^p)> bein S the P" th column of 
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Remark 2.5 The j.p.d. P (Z) in 1)2. 9|) is manifestly a symmetric function only of the 
eigenvalues of ZZ T , and thus, a symmetric function only of the singular values of Z. 

Remark 2.6 From the normalization condition (|2.24|) we obtain an alternative 
expression for the normalization constant (|2.23|) as 

1 _ 1 r dZ 

C ~ JdBf(tiBB T )\detB\ n ~ J [det(l + ZZ T )]( m + n )/ 2 ' ^' 2b > 

which is manifestly independent of the particular function f(u), in accordance with 
Lemma (|2.1|) . The integral over the matrix Z can be reduced to a multiple integral of the 
Selberg type [3J [7] over the singular values of the matrix Z, which can be carried out 
explicitly: 

dZ mn J\ r 

[det(i + zz T )]( m + n )/i = 71 2 W r(^) ' 

For particular choices of the function f(u), we can use (|2.25|) to derive explicit integration 
formulas. For example, the function 

/( M ) = 7rm (l+n)/2 ( 2 - 27 ) 

(i.e., the entries A{ a in (|2.1|) are i.i.d. according to a normal distribution of variance 1/2) 
satisfies 1)2. Thus, we obtain from p. 25)1 that 

2 n r ( m+j ) 

dBe -trBBT |detB|n = -Q \ 2 J (2 2g) 

i =1 ^(2/ 



Note that the integral on the left-hand side of 1)2. 28j) can also be reduced to a multiple 
integral of the Selberg type (this time, over the singular values of B), which can be 
carried out explicitly. The result is 



m£ -pr V 2 J 



Since this must coincide with ()2.28j) . we obtain the identity 



Example 1 For n = 1, i.e., the case where X and Z are m dimensional vectors, 1)27 
simplifies into the m dimensional Cauchy distribution 

= ( i + ztz)(^ • (2 - 30) 

This is so because for n = 1, the matrix ZZ T has m — 1 eigenvalues equal to 0, that 
correspond to the m — 1 dimensional subspace of vectors orthogonal to Z, and one 
eigenvalue equal to Z T Z. Thus, det(l + ZZ T ) = 1 + Z T Z . Eq. ()2~3U)) then follows by 
substituting this determinant into 1)2.22)1 . 
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3 The universal probability distribution of the product 
B l X of complex random matrices 

The results of the previous section are readily generalized to complex random matrices. 
One has only to count the number of independent real integration variables correctly. In 
what follows we will use the notations defined in the previous section (unless specified 
otherwise explicitly). Thus, consider a complex m x (m + n) random matrix A with entries 
Ai a (i = 1, • • • m; a = 1, . . . m + n). We take the j.p.d. for the m(m + n) entries of A as 

G{A) = /(tr^4t) = f(j2 \Aia\ 2 j , (3.1) 

with f(u) a given appropriate p.d.f.. From 6 

J G(A) dA = 1 (3.2) 
we see that f(u) is subjected to the normalization condition 

oo 

I u m(m+n)-l f ^ du = _J _ (3 3) 

« "2m(m+n) 

This implies that f(u) must be subjected to the asymptotic behavior 

u m(m+n) f ^ ^ q (3^ 

both as u — > and u — > oo. 

As in the previous section, we choose a partition 

A — > (B, X) (3.5) 

of the columns of A. Thus, tivL4 f = J2i.j \ B ij\ 2 + J2i, P \ x i P \ 2 = tiBB^ + trXX\ and thus 
(jHU reads 

G(B, X) = f(trBB^ + tiXX^). (3.6) 

We now define the random matrix Z = B~ 1 X. Our goal is to calculate the j.p.d. P(Z) for 
the ran entries of Z . The main result in this section is stated as follows: 

Theorem 3.1 The j.p.d. for the ran entries of the complex random matrix Z = B _1 X is 
independent of f{u) and is given by the universal function 

C 

= [det(l + ZZl)] m + n ' (3 " ?) 
where C is the normalization constant 



6 We use the Cartesian measure dA — d 2m ( m +") A — y[ dReAi a dlmAi a , with analogous definitions for 
dB and dX below. 
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Proof. The proof proceeds in a similar manner to the proof of Theorem (|2.1|) . The only 
important difference is that now 6(X) = UZi Up=i S(ReX ip ) 5(lmX ip ) = Up=i 5 {2m) {X p ), 
X p being the p-th column of X. One obtains 

P(Z) = J dBdX f(tTBB* + ti XX ] )5(Z - B~ l X) 

= J dB f[tvB(l + ZZ ] )B ] ]\detB\ 2n (3.8) 

where we have integrated over X. 

The m x m complex hermitean matrix 1 + ZZ^ can be diagonalized as 1 + ZZ* = UAIA^ 
where U is a unitary matrix, and A = Diag(Ai, . . . , A m ) is the corresponding diagonal 
form. Obviously, all Aj > 1, since ZZ^ is positive definite. Substituting this diagonal form 
into lUEEJ) we obtain 

P{Z) = J dB f{tiBUMA ] B ] )\AetB\ 2n = J dB /(tr BAB*) \detB\ 2n , (3.9) 

where we used the invariance of the determinant |det.B£/| = | det-B | and the invariance of 
the volume element d(BU) = dB under unitary transformations. 
As in the previous section we now rescale B as B = B\fK. Thus, 



det B = Vdet A det B , and dB = (det A) m dB . (3.10) 
Finally, substituting (|3.10|) in (|3.9j) we obtain 

dB ~ - t / |detS| x "" 



P(Z) = [ fjtrm 
1 ; J (detA)- M j V(detA)3, 



C C 



(3.11) 



(det A) m + n [det(l + ZZ^)] m + n ' 
where C is the normalization constant 

C = J dBf{tiBB^)\detB\ 2n (3.12) 

rendering 

' P(Z)dZ = 1 . (3.13) 



Finally, one can show, in a manner analogous to Lemma (|2.1|) . that C < oo and that it is 
independent of the particular function f(u). I 

There are obvious analogs to the remarks made in the previous section, which follow from 
Theorem ()3.1|) . which we will not write down explicitly. 

4 More on the distribution of solutions of systems of linear 
equations with random coefficients: extension of a result 
due to Girko 

The methods of the previous sections may be applied in studying the distribution of 
solutions of systems of linear equations with random coefficients. For concreteness, let us 
concentrate on real linear systems in real variables. 
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Consider a system of m real linear equations 

m+n 

^2 A iaia = h, i = l,...m (4.1) 

in the m + n real variables £ Q . With no loss of generality, we will treat the first m 
components of the vector £ as unknowns, and the remaining n components of £ as given 
parameters. Thus, we split 

«=(;) (4 - 2) 

where z is the vector of unknowns Zi = £j (i = 1, . . . m) and it is the vector of parameters 
u p = ip+m (p = l,...n). Similarly, we split the matrix of coefficients 

A = (B,X), 

where the m x m matrix B (with entries Bij) and the m x n matrix X (with entries Xi P ) 
were defined in (|2.6|) . Thus, we may rewrite (|4.1|) explicitly as a system for the Zf. 

Bz = b-Xu. (4.3) 

If we consider an ensemble of systems (|4.1|) . in which ^4 and 6 are drawn according to 
some probability law, the unknowns become random variables, which depend on the 
parameters u p . Girko proved jH] the following theorem for a particular family of such 
ensembles: 

Theorem 4.1 (Girko) 

If the random variables A{ a and b{ (i = 1, . . . m; a = 1, . . . m + n) are independent, 
identically distributed variables, having a stable distribution law with the characteristic 
function 

g(t;a,c) = e _c|t|a , 0<a<2;c>0, (4.4) 

then the random variables Z{ (i = 1, . . .m) are identically distributed with the probability 
density function 

oo 

p(C;a,{3) = | J rp(^—;ajp(r;a)dr, (4.5) 
o 

where p(r; a) is the probability density of the postulated stable distribution, and 



P=\i + 52K\ a j ■ ( 4 - 6 ) 

The ratios Zi/zj (i ^ j; i, j = 1, . . . m) have the density p{r; a, 1). 

In the special case a = 2 in (|4.4j) . the random variables Ai a and bi are normally 
distributed. For this case Girko obtained 
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Corollary 4.1 If, under the conditions of Theorem (|4.1|) . a = 2, then 



i.e., C follows a Cauchy distribution of width (3. 
When a = 2, the j.p.d. of the Ai a and the bi is 

G(A, 6) = (2n**) inil ^ ±11 e -5? ^ (48) 

which is a special case of the j.p.d.'s we have discussed in the previous sections. Thus, in 
the spirit of the discussion in the previous sections, we will study systems of linear 
equations (|4,1[) with random coefficients A and inhomogeneous terms b with j.p.d.'s of the 
form 

G(A, b) = f{t?A T A + b T b) , (4.9) 
with f(u) a given appropriate p.d.f. subjected to the normalization condition 

oo 

m(m + n+l) M 2, 

it* - 1 f(u)du = - . (4.10) 

^ ''m(m trt + 1) 

Our goal is to calculate the j.p.d. P(z;u) for the m unknowns Z\. We summarize our main 
result in this section as 

Theorem 4.2 If the random variables Ai a and bi (i = 1, . . .m; a = 1, . . . m + n) are 
distributed with a j.p.d. given by {4-^ , with f(u) being any appropriate probability density 
function subjected to \4- 10j ), then the random variables Zi (i = 1, . . . m) are distributed 
with the universal j.p.d. function 

(/3 2 + z T z)~ 

independently of the function f{u), where 



= VI + u T u (4.12) 
and C is a normalization constant given by 

C = J dA\detB\f{tvA T A). (4.13) 

Remark 4.1 Note that Theorem (|4.2jl generalizes the case a = 2 of Girko's result, 
Theorem ()4.1|) . from the particular f{u) ~ e~ u to a whole class of probability densities 
f(u), and moreover, it determines for this class of distributions the (universal) j.p.d. of 
the Zj's, and not only the distribution of a single component. Thus, it is an interesting 
question whether Girko's result could be generalized also to other ensembles of systems of 
linear equations as well. 
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Proof. By definition, from (|4.3|) . 

P(z;u) = JdAdbf(tiA T A + b T b)5(z-B~ 1 {b-Xu)) 

= J dAdbf(trA T A + b T b)\detB\5(Bz + Xu-b) 

= J dAf[tvA T A + (A£) T (AO]\detB\ , (4.14) 

where in the last step we integrated over the m dimensional vector b and used 
Bz + Xu = A£. 

The last expression in (|4,14|) is manifestly invariant under 0{m) x 0(n) orthogonal 
transformations 

P(O lZ ;0 2 u) = P(z;u), O x G 0{m) ,0 2 G 0(n), (4.15) 

due to the invariance of the measure dA = dB dX = d{BO\) d{X0 2 ), the invariance of the 
determinant |detS| = |det-B0i|, and the invariance of the trace 

tr A T A = tr (BOxfiBOi) + tr (X0 2 ) T {X0 2 ). With this symmetry at our disposal, we 
may simplify the calculation of P(z; u) by rotating the vectors z and u into fixed 
convenient directions, e.g., into the directions in which only z\ and u\ do not vanish: 



with 



zP=z 5 n , u^=u 5 pl . (4.16) 

T„\l „, _ f„.T„ - 1 



z = (z 1 z)2, u = {u 1 u)^. (4.17) 
Thus, we obtain 

P{z-u) = J dAf[trA T A+(Bz^ +Xu^) T {Bz^ + Xu^)]\detB\ 

= J dBdX f(S)\detB\, (4.18) 

where 

m n 

S = E B J B J + E X p X p + 0- + z o) B f B i + 2u o*o Bf X l + (1 + ug) XfXx , (4.19) 

3=2 P=2 

in which Bi is the i-th column of B, and X p is the p-th column of X. 
The bilinear form involving B\ and Xi in (|4.19f) may be diagonalized as 

(l + ^)( a0B ^ 0Xl ) 2 + ( "° B ^| X ' ) 2 . (4.20) 

where we have used + u o = We now perform a rotation in the I?i — X\ plane, 
followed by a scale transformation of the first term in Q4.2UJI , thus defining 

B[ = (l + e& Z0Bl t^ X \ X{ = UoBl Z^ Xl , (4.21) 
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such that d m B[ d m X[ = (1 + £ T £)^ d m B 1 d m X 1 . We will also need the inverse 
transformation for B\ 

BliB'^Xi) = * ( ZoB j +u X[) (4.22) 

in order to express the matrix B in terms of the primed column vectors: 

B = {B 1 (B' 1 ,X[),B 2 ,...,B m ) . (4.23) 

Thus, using dOH]) - and the trivial fact that dA = ]T™ i d m B t U7= P d m X p , we 

obtain 

«) = / dA /(tr ^ } |d6tjg| ' (4 - 24) 

where we have removed the primes from the integration variables. We are not done yet, 
since B\, the first column of B, depends on zq and uq. To rectify this problem, we note 
from (|4~22")) that 



BiiB^Xt) = J 1 (B 1 cos + Xt sin ) , cos fl = ^^^=. (4.25) 

V 1 + ^ */^(i + u§) 



Thus, performing one final rotation by an angle 6 in the B\ — X\ plane, which leaves, of 
course, dA and trA T A invariant, we see that in terms of the rotated columns 

/ l+ u 2 

|det.B| = y |detS|, and thus, finally, we obtain that 



^(*;«) = ; V - / dAf(tvA T A)\detB\ , (4.26) 

which coincides with (|4~TTj) . due to (|4~T2"|) and (JUTS}. ■ 

Remark 4.2 The fact that the integral C = J dA |det£>| f(trA T A) is convergent and 
independent of the function f(u) can be proved by decomposing A into its singular values 
[SI HI , essentially in a manner similar to our proof of Lemma (|2.1|) , but with slight 
modifications in (j2.13j) and (|2.14j) due to the fact that A is a rectangular matrix rather 
than a square matrix. We shall not get into these technicalities here, which the reader may 
find in Note, however, that C may be determined from the normalization of P(z; u): 



1 



d m zP(z;u) =1 = CP [ m+1 = ■ 



Thus, 



r(^±i) 



C = = ^ • (4.27) 

7T 2 "Jm+1 
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Remark 4.3 We note that for n = (i.e., when A = B) and u = 0, Q4.3J1 degenerates into 
Bz = b, which is precisely the case n = 1 (and z = Z) m. the conditions for Theorem (|2.1|) . 
which we analyzed in Example Thus, (|4.11|) . evaluated at n = 1 and u = must 
coincide with (|2.3(Jj) . as one can easily check it does. 

Since Theorem (|4.2j) states the explicit form (|4.11j) of P(z;u), we can now use it to derive, 
e.g., the probability density of the distribution of a single component z% and that of the 
ratio of two different components, mentioned in Girko's Theorem (|4.1jl : 

Corollary 4.2 The m components z% are identically distributed, with the probability 
density of any one of the components Zi = C, given by Q4.7|) of Corollary ()4.1|) . 

Proof. That the Zj are identically distributed is an immediate consequence of the 
rotational invariance of P(z; u) in ()4.11|) . The proof is completed by performing the 
necessary integrals: 



Thus, from (|4.27|) we obtain the desired result that p((; (3) = ■ This result, should 

have been anticipated, since the universal formula Q4.11[) holds, in particular, for the the 



Finally, we have 

Corollary 4.3 The ratios Zi/zj (i ^ j; i, j = 1, . . . m) have the density P(r) = p(r; 1). 

Proof. The ratio Zi/zj is dimensionless, and thus its distribution cannot depend on the 
width f3, which is the only dimensionful quantity in l|4.11|) . The proof amounts to 
performing the necessary integrals, e.g., for the random variable z\j z<i'- 





(4.28) 



Gaussian distribution (|4.8jl . 



P{r) = f d m zP{z;u)s(r-—) 




z 2 \dz2 ...dz 

in 



Z2J 



m+l 

2 




(4.29) 
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